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ELLIPTICAL ORBIT-PERFORMANCE 
COMPUTER PROGRAM 


by T. R. MYLER 
VOUGHT CORPORATION 


SUMMARY 


This report describes and presents a FORTRAN coded computer program 
which generates and plots elliptical orbit performance capability of space 
boosters for presentation purposes. The program requires input data from a 
trajectory simulation which defines the booster’s velocity capability as a 
function of insertion altitude and payload weight. The Elliptical Orbit 
Performance computer program manipulates the velocity-altitude-payload weight 
data to obtain apogee altitude-perigee altitude-payload weight data and 
generates a computer plot. Included in this report are program theory, user 
instructions, output definitions, subroutine descriptions and detailed 
FORTRAN coding information. 
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1.0 INTRODUCTION 


A common method of presentation of orbital performance capability of 
space boosters is to show apogee and perigee altitude as a function of payload 
weight. Typically, apogee and perigee altitude data are calculated from 
parametric data of altitude and velocity at orbit insertion. The booster’s 
velocity capability as a function of altitude and payload weight at orbit 
insertion are commonly calculated by a computer program which simulates the 
booster flight. Thus, based upon the parametric results of a trajectory 
program which simulates a specific space booster, the orbit insertion data can 
be manipulated to produce parametric data of apogee and perigee altitude as a 
function of payload weight. The mechanization of this manipulating process 
resulted in computer program Elliptical Orbit Performance (acronym ELOPE) and 
is described herein. 



2.0 DEFINITIONS 


Flowchart conventions used in this report are as follows 



Process 



Input/Output 


Subroutine 



Decision 


Subroutine Call 



5.0 PROGRAM DESCRIPTION 


This section describes program theory, input instructions and output 
definitions. 


3.1 Program Theory 


The purpose of computer program ELOPE is to calculate parametric data 
in apogee altitude, perigee altitude and payload weight and generate a 
computer plot. The technique used in ELOPE to obtain this plot includes 
interpolation of data and the solution of a two body energy equation. 


Parametric data in altitude, velocity and payload weight at orbit 
insertion must be input to ELOPE. These data can be obtained by use of a 
trajectory program, such as NEMAR of Reference (l), which simulates booster 
flight. By calculating trajectories with various values of payload weight and 
vehicle pitchover rate, the data map depicted below can be obtained. 


Insertion 

Altitude 

(Perigee 

Altitude) 



The above information is interpolated at the perigee altitudes of 
interest for each payload weight. Perigee is assumed to be at the point of 
insertion since insertion at perigee results in maximum vehicle performance 
and since orbiting vehicles commonly insert at perigee. The interpolated 
velocity-payload weight at each specific perigee is itself interpolated for the 
range of payloads weights defined by input. The resulting values of perigee 
velocities are converted to apogee altitudes by the following relationship: 


h 


a 


2 

2 - Vp^ 

rp GM 


- Tp - re 


where 


= perigee radius = hp + re 
Vp ~ perigee velocity 


= Earth radius 

GM = Earth’s gravitational eonstant 
= 1.4076576 X 10^‘’ft^/sec^ 
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The parametric orbital performance plot is defined by repeating the 
above described process for the range of perigee altitudes defined by input. 

Following the payload weight calculation cycle for a specific perigee 
altitude, payload weights are calculated which correspond to circular and 
Ehrth escape orbits. Circular orbit velocity is obtained from: 


Vc = 


Earth escape velocity is obtained from: 


Ve = 


The previously defined velocity-payload weight data are interpolated at 
the above two velocities to define the circular and escape payload weights. 

When the payload weight is desired for a specific orbit, in lieu of 
parametric performance, the velocity-payload weight data is calculated for the 
perigee altitude of interest as previously described. The velocity required 
for the specific orbit is then calculated from: 




Vf = GM ( 2 - 1 ) 

fp 3. 


where a = semi-major axis = + Tp 

2 


r^ = apogee radius = h^ Tq 


This perigee velocity is used to interpolate the velocity-payload 
weight data for the payload weight corresponding to the specific orbit. 


3«2 User Instructions 


ELOPE uses a modified FORTRAN NAMELIST for inputting data which 
provides the user with readability and simplicity of use. 

The following rules apply to NAMELIST used by ELOPE: 

1. First card of a data group or case is ^INPUTD beginning in coliimn 
2. Blanks are not allowed. 

2. Last card of a data group or case is $END beginning in column 2. 
Blanks are not allowed. 
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Blanks may not be used within names but may be used elsewhere. 

4* Variable names are followed by an equal sign which is followed by a 
value which is followed by a comma, e.g., WEIGHT = 323. 07? 

5» Only columns 2-72, inclusive, are used. 

6. Titling information may be input by the appropriate title names, 
e.g. , TITLE1= ELLIPTICAL ORBIT PERFORMNCE - VAFB LAUNCH 
TITLEl must begin in column 2. 

7* Any number of names and values may be on a single card or line. 

8. Complete data arrays are input in the following form: 
name = value , value , value , . . . , 

Data values may be continued on the next line, but the last 
character on every line must be a comma, excluding title cards. 

9. Repeated data values may be input by using a repetition factor and 
an asterisk, e.g., VI = 38050, 36525, 2*31510, 28450, 

10. One or more specific elements of an array may be input, e.g., 
WEIGHT(3) = 200, 300, 

Subsequent data cases are allowed by providing additional sets of 
NAMELIST data. All input data is retained for subsequent cases but can be 
changed by inputting new values. 

A sample data case is included as Appendix A to exemplify data case 

setup . 


Execution of ELOPE requires that the CalComp pen plotting facility be 
available to the computer job at the time of program load. This facility 
consists of the CalComp 763 pen plotter hardware and the CalComp Basic 
Software Package, Reference (2). ELOPE generates three plots per case if 
selected by input. A non-zero value of IPLOT results in plotting on graph 
paper which has a perforation size of 11 x 17 inches, a gnd size of 9 7/8 x 
15 inches and gnd tjrpe of 10 divisions per centimeter. At the Vought 
installation, this paper is identified as CAL32. A non-zero value of LPLOT 
results in plotting on 4 cycle semi-log graph paper which has a perforation 
size of 11 X 8 l/2, a grid size of 10 x 7 and gnd type of 2 l/2 inches per 
cycle and 10 divisions per inch. At the Vought installation, this paper is 
identified as CAL44* Additionally, this same data is plotted on no-gnd paper 
and may be scaled down in size for use in vugraphs. This plot is placed on 
paper identified as CAL36 at the Vought installation. 

Definitions of specific NAMELIST inputs to ELOPE are shown below. 
Default values are shown when they are set by the program prior to reading 
input data. 
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NAMELIST Input Definitions 


ALIMAX 

Maximum perigee altitude used for parametric 
output. Units according to lOPT. Use when 
lOPT = 1 or 2. 


ALIMIN 

Minimum perigee altitude used for parametric 
output. Units according to lOPT. Use when 
lOPT = 1 or 2. 


APOGEE 

Apogee altitude of single orbit case. Units 
according to lOPT. Use when lOPT = 3 or 4. 


DEL ALT 

Increment in perigee altitude for parametric 
output. Units according to lOPT. Use when 
lOPT = 1 or 2. Maximum number of altitude 
points is 50. 


DELWT 

Increment in payload weight for parametric 
output. Units according to lOPT. Use when 
lOPT = 1 or 2. 


FACT 

Ratio of plot size to normal plot size. When 
greater than zero but less than or equal to 
one, the semi-log graph is plotted on no-grid 
paper. This plot resides on local file name 
PLT3. (0. built-in) 


lOPT 

Output data option 

= 1 Input ALTMIN, ALTMAX, DELALT, WTMIN, 

WTMAX, DELWT and output parametric data 
in n.mi. and lbs. (1 built-in) 

= 2 Input ALTMIN, ALTMAX, DELALT, WTMIN, 

WTMAX, DELWT and output parametric data 
in km and kg. 

= 3 Input APOGEE, PERIGE and output single 
orbit only in n.mi. and lbs. 

= 4 Input APOGEE, PERIGE and output single 
orbit only in km and kg. 

- 

IPLOT 

Non-zero value produces a CalComp plot of 
altitude as a function of velocity. Local 
file name of plot is PLT2. (0 built-in) 


IPRNT 

Output frequency control. Parametric data is 
calculated at DELWT intervals from WTMIN and 
printed at the IPRNT frequency. (1 built-in) 


IRAD 

Input data option 

= 1 Input R1 - R15 as radius in feet. 

(1 built-in) 

= 2 Input R1 - R15 as radius in n.mi. 

= 3 Input R1 - R15 as altitude in n.nm. 


KASE 

Case number. 
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LPLOT 


PERIGE 


PLABELl - 
PLABEL9 

PTITLEl - 
PTITLE4 

HEARTH 


R1 - R15 


TITLBl 
TITLE2 
VI - V15 
WEIGHT 

WTMAX 

WTMIN 

XINC 


Non-zero value produces a CalComp serai-log 
plot of apogee altitude as a function of 
perigee altitude and payload weight. Local 
file name of plot is PLOT. (O built-in) 

Perigee altitude of single orbit case. Units 
according to lOPT. Use when lOPT = 3 or 

labels placed on upper right side of apogee- 
perigee plot. Maximum of 30 characters each. 

Titles placed at top center of apogee-perigee 
plot. Maximum of 40 characters each. 

Ehrth radius used to calculate altitudes, ft. 
(20925741. built-in) 

Tables of radius or altitude (according to 
IRAD) for each WEIGHT. Input in increasing 
order. Minimum of 4 and maximum of 10 values 
per table. Enter 0. after last value of each 
table if less than 10 values are input. 

Minimum of 4 tables. 

Title printed at top of each page. Maximum of 
72 characters. 

Title printed at top of each page. Maximum of 
72 characters. 

Tables of inertial velocity in fps for each R1 
- R15- Maximum of 10 values per table. 

Table of payload weights in lbs. 

Corresponding tables of radius and velocity 
must be input for each WEIGHT. Maximum of 15 
values. 

Maximum payload weight used for parametric 
output. Units according to lOPT. Input when 
lOPT = 1 or 2. 

Minimum payload weight used for parametric 
output. Units according to lOPT. Input when 
lOPT = 1 or 2. 

Increment value of x axis major divisions of 
apogee-perigee plot. Units according to 
lOPT. There are 7 major divisions on the x 
axis. (lOO. built -in) 
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Output Description 


The NAMELIST input data is listed verbatim as read. This list provides 
a quick check of the input data for format correctness and validity. 
Additionally, the parametric data of weight, altitude and velocity are output 
in a different format than input for inspection purposes. 

Subsequent pages provide parametric data of perigee altitude, apogee 
altitude and payload weight. For each perigee altitude, as defined by the 
input, apogee altitude is calculated at each payload weight increment from the 
minimum value to the maximum value. The resulting values, in both English and 
metric units, are output according to the value of the input IPRNT. Following 
the parametric data of each perigee, the payload weight corresponding to a 
circular orbit and to escape velocity are shown. A sample data case is 
included as Appendix A. The plots resulting from this data case are also 
included in Appendix A. 

For single orbit cases, data for the specific orbit are shown in lieu 
of the payload-apogee altitude parametric data. 
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4.0 SUBROUTINE DESCRIPTIONS 


This section provides a brief explanation of each subroutine of ELOPE. 


4.1 ELOPE (Main Program) 


The main program initializes the input data defaults; converts input 
data to internal units; calculates parametric data of apogee altitude, perigee 
altitude and payload weight; calculates payload weight at circular and escape 
velocities; writes parametric data on disk units 1 and 2 for subsequent 
plotting; calculates single orbit payload weights; and outputs the results. 


4.2 INPUT 


Subroutine INPUT reads input data in a modified NAMELIST format. 

Titling information on title cards are placed in appropriate arrays for use by 
the main program. Non-title cards are written on disk unit 8 for a FORTRAN 
NAMELIST read from the main program. 


4.3 INTER 


Subroutine INTER is a second-order interpolator of two variables. It 
selects the four closest data points to the desired value of the independent 
variable and interpolates or extrapolates for the value of the dependent 
variable. 


4.4 PLOTLG 


Subroutine PLOTLG produces a CalComp semi-log plot of apogee altitude 
as a function of perigee altitude and payload weight. The ordinate is fixed 
at a four cycle logarithm scale from 100 to 1,000,000 n.mi. or kilometers. 

The abcissa scale is determined by input data but is limited to seven major 
divisions in length. Plotted data is taken from disk units 1 and 2, which are 
written by the main program. The data plotted is scaled according to an input 
multiplier in order to decrease the physical size of the plot. This 
multiplier is one when plotting on grid paper. 
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4.5 PLOTLR 


Subroutine PLOTLR produces a CalComp plot of insertion (perigee) 
altitude as a function of insertion velocity. These data are input values and 
are plotted for inspection purposes only. The ordinate is a fixed scale from 
0 to 1200 nautical miles. The abcissa is a fixed scale from 18000 to 37000 
feet per second. Also plotted on this graph are the circular and escape 
velocity lines. 
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5.0 PROGRAM CODING 


This section presents details about the program coding. Included are 
flowcharts of each subroutine, FORTRAN listings of each subroutine and 
definitions of the FORTRAN variables. The information presented in this 
section is intended to be helpful in developing a thorough understanding of 
ELOPE and in making modifications to the program. 


5.1 Subroutine Flowcharts 


Flowcharts are presented in Figures 5.1 through 5.5. 
conventions used in these figiires are defined in Section 2.0 


Flowchart 
of this report. 


5.2 FORTRAN Listings 


ELOPE is coded in FORTRAN IV, Reference (3), on the CDC CYBER 175 
ccanputer with the NOS/BE 1.4 operating system. Listings of the FORTRAN coding 
are presented in Appendix B. 


5.3 FORTRAN Variable Definition 


Definitions of the FORTRAN variables are shown below. This information 
is usually used only when making modifications to the program. 


Variable 

Definition 

A 

Semi-major axis, n.mi. 

ALTMAX 

Input value 

ALTMIN 

Input value 

APOGEE 

Input value 

DELALT 

Input value 

DELWT 

Input valve 

DH 

Perigee altitude increment 

DW 

Payload weight increment, ; 

FACT 

Input value 
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Variable 


Definition 


FTNM 

GM 

H 

HA 

HAMAX 

HAMET 

HF 

HH 

HHMET 

HMAX 

HMIN 

HNEW 

HP 

HPMET 

lOPT 

IPAGE 

IPLOT 

IPRNT 

IRAD 

KASE 

LABEL 

LPLOT 


Feet per n.mi., 6076.11549 

Earth's gravitational constant, 1.4076576 
X 10l6 ft3/sec2 

Table of altitudes from R1 - R15i n.mi. 

Apogee altitude of single orbit case, n.mi. 

Maximum apogee altitude plotted - 
1,000,000 - n.mi. or km 

Apogee altitude, km 

Perigee altitude, n.mi. 

Table of perigee altitudes of parametric 
output, n.mi. 

CXirrent perigee altitude, km 

Maximum perigee altitude of parametric 
output, n.mi. 

Minimum perigee altitude of parametric 
output, n.mi. 

Current perigee altitude, n.mi. 

Perigee altitude of single orbit case, 
n.mi. 

Perigee altitude, km 

Input value 

Page number 

Input value 

Input value 

Input value 

Case number 

Qiaracter data of information input via 
LABELl - LABEL9 

Input value 
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Variable 


Definition 


NC 

NDIM 

NH 

NNW 

NREC 

NWT 

PERIGE 

PL 

PLM 

R 

REAR 

RE 

REARTH 

RF 

R1 - R15 
TITLE 

TITLEl 

TITLE2 

V 

VEL 

VELC 

VELE 


Tables of number of non-zero values in the 
R1 - R15 tables 

Number of permissable values in R1 - R15 
and VI - V15 tables, set to 10 

Number of perigee altitudes of parametric 
output 

Number of payload weights of parametric 
output 

Number of records of data written on disk 
unit 1 

Number of non-zero values of the WEIGHT 
table 

Input value 

Current payload weight, lb 
Current payload weight, kg 
Array of R1 - R15» input units 
Circular orbit radius, ft 
Earth radius, n.mi. 

Input value 
Perigee radius, ft 
Input values 

Character data of information input via 
PTITLEl - PTITLEil 

Input value 

Input value 

Array of VI - V15, ft/sec 
Current velocity, ft/sec 
Circular orbit velocity, ft /sec 
Escape orbit velocity, ft/sec 
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Variable 

Definition 

VI 

Table of velocities at the current perigee 
altitude, ft/sec 

VMAX 

Velocity corresponding to an apogee of 
1,000,000 - ft/seo 

VV 

Tables of velocities for each perigee 
altitude of paremetric output, ft /sec 

VI - VIO 

Input values 

WEIGHT 

Input value 

WI 

Table of payload weights at the current 
perigee altitude, lb 

WMAX 

Maximum payload weight of parametric 
output, lb 

WMIN 

Minimum payload weight of parametric 
output, lb 

WNEW 

Current payload weight of parametric 
output, lb 

WTC 

Circular orbit payload weight, lb 

WTCMET 

Circular orbit payload weight, kg 

WTE 

Escape orbit payload weight, lb 

WTEMET 

Escape orbit payload weight, kg 

WTMAX 

Input value 

WTMIN 

Input value 

WW 

Tables of payload weights for each perigee 
altitude of parametric output, lb 

XINC 

Input value 

XKG 

Kilograms per lb, 0.45359 

XKM 

Kilometers per n.mi., 1.852 
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Figure 5.1 

FLOWCHART OF MAIN PROGRAM ELOPE 










Figure 5.2 

FLOWCHART OF SUBROUTINE INPUT 
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Figure 5.3 

FLOMCHART OF SUBROUTINE INTER 
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Figure 5.4 

FLOWCHART OF SUBROUTINE PLOTLG 
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Figure 5.5 

FLOWCHART OF SUBROUTINE PLOTLR 
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409111''^ 

'•TICiL '*'5'’JT 3C 



tNC-LIS:-* 1 ' 

"’ITS 


0 ^'■93 

OLHor 






r 

AS 2 i 

PAG? 1 


! '■' E - ALT. 

- 1 

.'u.C 

* 1 ** 

5.2 •< ■'1 



IME-i it:tv 

TATA (Oi- 

EC^ 

FO? ‘-CCU'A 

C^) 




,- 5 V L ^ i I '■ 

J T I ?- 







-T. f 1 as V 

= L. »P = S 
















a*---';.- 







.t 

4 , 3 ; . e 







3 ” . . L 

1 







3 C'- . •: 

3<':i4>^ 







4 : : . ; 

?P?®0 .1 







K ' 

S6P42 ,q 







‘•OC.C 

2f *'17.6 







ogp rr^zc 

payload 


APOGEE 


p E I G E E 

PAYLOAD 

APOGEE 

ALT..N‘.MI. 

w T . » L 3 S 


A L T . , . M I 

. 

VEL,»FOS 

■W T . » K C- 

ALT.jKM 

1 U • C 

120.0 


67097.0 


35262. B 

54,4 

124263.6 

100 . c 

140.0 


323B1.5 


344 = 2 . 5 

63.5 

59970.6 

1 50 • 

IM.O 


20356.2 


3377« .9 

72.6 

33625.7 

l'‘0.0 

ISu.O 


15056. e 


33124,6 

31.6 

27385,2 

100.0 

200. C 


116C9.1 


32521.1 

= C.7 

21500.0 

lOC.C 

220.0 


932C.1 


31 = 6 = . 1 

= = .6 

j.7260.8 

1 jC . c 

240 .0 


7546.5 


31477.1 

136.0 

14161, 3 

1 rc.c 

2 4^.G 


*'"9 c 


31 > ) 2 . 4 

117.0 

1 1-3 3 9. 7 

i ‘ ^/ 9 



5 4 1C 


2 .560 .= 

1 2 7 . C 

r:L36.1 

I. 'w • L' 

2 •■> 


4o3<4.7 


3 014 3 .6 

136.1 

-*5 = 2.7 

1 "'*.0 

3 20.0 


399^.9 


2 = 755.4 

145.1 

7391.1 

1 1 ‘ • c 

3 4 J ^ 


5 4 to .7 


2 = 344.7 

154.2 

63 f=7.0 

lOC.C 

3^C. C 


-) Q j -j , 


24034.2 

163.3 

5534.9 

: ; : . *: 

•3 J 

^ > • -« 


PBc-a. f 


2 ' 70 i , f 

17 2.4 

4 c r 1 . 0 

.V • i 

4 .3. w 


2247.5 


2 4 -5 = . 1 

13 1.4 

4162.3 

IC‘C.0 

42J. 0 


lo3»- . <. 


.4 

19 3.5 

3 5 « = . 7 

3 j u ■ w 

4 4 ^ , j 


1 ^ “D V . ^ 


2 7773 .= 

x = 4.4 

3075.3 

^ • i. 

-.60 .C 


1414.4 


Z74=4^n 

2 ■: 4 . 7 

741 = . 5 

1 W V* * Vf* 

4^3.: 


il<=3 . 7 


: 7 2 3 ? . 1 

717.7 

2 2 A 0 .4 

i*:c.c 

■jCO.C 


494.0 


r-=6r .= 

226 .r 

1442.6 

1 L C . 0 

= 20.c 


«lr .1 


26 713 . 7 

235 .= 

15 09. 7 

1 "J.C 

540.0 


63 J. • 4 


2c;47'^ . 1 

744.0 

1207,3 

2 V v« • C 

560.0 


5^3.1 


2624': .: 

754 . 1 , 

= 3x .7 

100. c 

5«C.0 


3^4.0 


P4C27. 6 

7 = 3,1 

679,6 

^ s ^ 

^ J . Iw 

;.o 


243.1, 


17.-4 

773 .2 

4 4 c , 2 

1 V, ' # • u 

^2 0. C 


127.0 


25614 . 3 

231.2 

235.1 

C!-CIjLAP ’AYL'"A0= ol-i."' 

LBS 

= 2 ".3. 3 i<r- 


CT=:iJL44 V:L. 

= 26507,? 

Foe 

•ES'aPE PftVLQAr 

= 100.6 

L-5 

= 4 -. 7 i/r 


= S''iP- \'^L. 

= 9 

c ? e 


- A4 


f'r J) in iJ uj sO o H* -g to oo 



LLI = TICiL r-ilT 


tMr,LlS-J LNITS 


C^Kisns AT If'N 


>' i.\r 


p f fiC 


P 1 : T r P P A L T « = : 


TEROL AT I' 


"'.“T, X 37^,4 ri ‘-\ 

>Ti (CHEC*^ '^"’5 ACrijAACY) 


YL ■ AC 

IV j,-:rTinv. 

. . 1 3 s 

V t E . , ? P S 

-» • 

p - P 7 . . 


? -' ; p . 

^ C " 

^ A ^ 

^ ^ J ^ 

. ^ ^ • o 

ZZ* • ( 

311..S.7 

3C"*c 

2 c? -3 2 , 4 

4^*0, 

27^3j.2 

. L V • V. 

2 f 3-?^ , 4 

^ Vj V* • U 

2 4^ 1 5 . 5 


PERIGEE 

PAYLCAO 


AP3GE= 

PERIGEE 

PAYLOAD 

APOGEE 

ALT.f 

/I T ..LBS 


ALT .^^■.’^I. 

VEL.i EPS 

W T . > K G 

ALT. >KM 

200 . C 

ICD.O 


2C4'3PB ,4 

35351. C 

48,4 

370457 . C 

2:o.c 

12 0. C 


47260.1. 

34442,4 

56.4 

87536.7 

2 00. r 

140.0 


28L13.5 

33665.0 

63.5 

68362.2 

200. C 

IfO.O 


17M4.6 

32923 .fc 

72.6 

32622,5 

2 00. C 

180.0 . 


12991.7 

32261.3 

81.6 

26C60.7 

2 o 0 .0 

200.0 


i'j 122.1 

31659.3 

90.7 

1=»766.2 

2CC.G 

2 20.O 


=> 164... 

Slice. 7 

99.3 

15119.8 

200. P 

240. C 


8"'.-6.1 

30587.0 

i C 8 . 9 

12417.8 

2 '. ? . C 

2 6g. 0 


EBCc.6 

3vU4.l 

1L7.9 - 

- 10344,0 

2 ' ( v"' . 0 

28 0.0 


472t.T 

2<^653. 7 

127 .C 

9752.6 

2 .C 

3CJ.0 


4^23. t 

2^232.6 

136.1 

7651.5 

2 00.0 

320.0 


3434.4 

2B33C. 8 

14 5,1 

53 oO . 6 


3 4 : . 0 


? G3 3 , ^■ 

2=451.2 

154.2 

6644. ■» 

i , r 

3-0.0 


Cfl’ .2 

2=:^3 .7 

1 c p • 3 

4443.7 

C_ V 

3 ? 0 . J 


21-4,7 

27753.0 

172. 4 

3990.5 

^ JO • c 

4 w 0 .0 


1337,6 

2'’43l .2 

i"'1.4 

34c? .2 


420.0 


^554. i. 

27118 .0 

l-PC .5 

287?. 1 

z^c.c 

440.0 


13-'2.3 

2 85 1 c , 7 

199,4 

2611.9 

Z jC 

480. C 


10^7.? 

2 6 = 3 4 , r 

2 !- P . 7 

1 = 96.1 

2 :‘". 

48 0. C 


P74.C 

26261,6 

217,7 

1420.4 

200 *C 

5CO.O 


691.9 

2599 c. A 

226 .8 

L 2 S 1 . 4 

2-0 #C 

5 20.C 


f 25 . 4 

25747.5 

23*= .9 

972.1 

2 y 0 • C 

?40 .0 


373.2 

2 5 5 C 5 . 5 

2 44,9 

491.1 

2 V? • C 

5 6 -J.C 


233.2 

25271.6 

254. c 

631,9 

CIS^CULA=^ '>AYLr^ 

AC= 58h*9 

LBS* 

256.2 4C 

CIRCLLA? VEL 

.* 25214.5 

EPS 

o^yl:?ad 

“ R3 . c 

LBS* 

4 2 . <r- 

ESCAPE ''EL. 

* ;5-35=.7 

c 0 5 


- A5 



G'‘'l'T 1 'U’lr- '/ = -TrL = 


vniJG^T C?i5p?9ATIC‘'J 

LLl'=T''CiL 33‘«iT PrppO? 

rNSLI--* L'-ITS 

^lCF= 

CiSE 1 -ir- 3 


T M T r. c = 

= alt.* 3u:.> * 555. c 

'l.iTITM r^TA (GhE''k -C' iGGM^ACY) 

= M Y L ■> i 0 

: 4jcoT:<-f 


VEL. >Fos 

* f 

“ * • 

3 - 

• . 

3 5:. .5.0 

1 = ' . : 

33512.2 

. • f • V 

3 .211.W 

3'^', * !. 

2P’C5.^ 


2645.3. C 

50-uC 

25C12.4 

fcCu.o 

238J3.9 


D.-C JCP = 

0 tYL3A3 

iFGGEc 

P£9 IGEE 

PAYLOAD 

APOGEE 

4LT.,N.f"I. 

WT., LBS 

ALI.t^.-^I . 

VEL. ,FPS 

•n . #K6 

alt .,km 

ao'-.c 

100. u 

91042 . 1 

34502.4 

45.4 

158621.1 

3 ‘3 . 0 

120.0 

35iC 7.C 

3356c .2 

54.4 

55019.3 

3 o o . 

19G.0 

21235.2 

327cC.6 

52.5 

39329,3 

300.0 

150.0 

14644.5 

32051.2 

72.5 

27492.2 

300. C 

leo.c 

3113c. 5 

3138C. 3 

61.5 

20524. g 

3CC.C 

200. U 

6745.1 

3077t .0 

9C.7 

16197.9 

300.0 

220.0 

7374.3 

30211 . 0 

99.8 

13101.5 

3 : : . 

240. J 

53 J7.5 

295P2.2 

iOR .G 

1C755.5 

?rc.c 

?60.0 

4.^32 .2 

2=^19 0 .G 

1x7. G 

?G49 . 3 

^ ' r\ 

^ J\, 

2BC. C 

40 55.3 

2 7 ? 2 . c 

177. C 

7515.0 

3 V. 0 . C 

w w Lrf • Stf 

342S .5 

2 6 2 C 4 . 5 

13c .1 

634Q .q 

300.0 

3 20.0 

2R47. r. 

. : 

1 4 R . 1 

53c5.2 

3^:.: 

3 4'j.C 

2 4 46 .1 

2 7 5 : ? . 7 

-54.2 

45 3 3 .7 


3-0. C 

2053 . - 

27143.3 

1 r 3 . ? 

3^22.2 

3*''C.^ 

35 0.0 

1721.2 

2c7G7.i 

172.4 

32v.5.2 

300. C 

4 jo. : 

' 44f. 1 

7 fjtt f ^ j 

i^l.^ 

2567.1 

3 ^r.c 

42J.C 

1L.--J .C 

2 5 1 0 . 7 

192.5 

2 1 b r . 3 

3 ’' v.O 

44.^. 0 

^46.3 

2 - 6 4 7 . c 

1 G 9 , '■ 

1756.3 

2 1 C . L 

453. C 

74 j. 7 

2 5 5: 7.5 

2:.R .7 

1371 .4 

3CC.C 


553.5 

25279.5 

217.7 

1^25. i 

3 . 0 

5V..3.L 

333.5 

25012.4 

2 2 0 . R 

71 C. 4 


flRCULAP 3AYLQAC* 51G.5 L8S» 131. t KG 


p^rftpf: OiYL"fin 


U ^ 
X y* 


CI'^C'JLAR V-:l.» 2^-7=,: 
2 SGAC’? VrL. = 


■ 5 $ 


L?S = 


3 ^ . 4 


3 51 2 ?PS 



!.ii''!C 4 vjHicis 

= LLi = t:C'.L 3c?cnp ^ Ef'GLiJ-^ '.'^ITS 


o ^LT.» '-0C.U 4.-1. = 


p,XE-®'^L ttT’'?^ •■'ATA (CHECK F^5 ACCJ'^ACV) 

- A Y L “ A r !“" J ‘^C T I 

'^T. f 1*^5 *V = L*. 

? 7 iC^.O 

- M 1 -i* 0 , ^ 

r •- ? 7 , c 

p , Q 

? 73 a ?, 1 

2 'i.Qi.e 
23 cc ^ ,5 
227 ’^ 0 , u 


0E5IGEE 

payload 

— 

APQGEE 

PERIGEE 

payload 

APOGEE 

ALT.^M.'*!. 

’*T., LBS 


A L T . 1 N , ^ I , 

vel.,fos 

iJ T . » K G 

ALT., KM 

'maO.C 

IvO.O 


55H65.3 

3363S.8 

4= .4 

10275P.8 

ACC.C 

120.0 


26^<?1. 6 

32714.2 

84.4 

49998.2 

A > 0 . b 

1 HO . 0 


17339.2 

319CC..2 

0 3 . 6 

322C4.9 

AuC.C 

160.0 


12489. t 

31162.3 

72.6 

23127.7 

'*00.0 

l = C.C 


94b3.p 

304F2.P 

81.6 

17873.2 

AuC.O 

2 \/ 0 . 0 


749C.7 

2<5864.0 

90.7 

13372.0 

A0).O 

22*). 0 


6061. c 

29296.9 

9=.F 

11226.7 

aCC.C 

24..0 


4963.1 

2 3 7t. J.7 

1<,8 .9 

9191.7 

A T , 0 

c '• 0 . c 


AlCr-.T 

28262.1 

117.9 

7608.6 

Ai.-: .0 

23'J.f 


3-2'.. 4 

27797.3 

127.0 

''334.6 

‘ 0.0 

3oC .Cl 


23=7.6 

27362.1 

13o .3. 

=292.2 

A V, V .0 

3 2*0.C 


2 3 7-3.2 

269i»5 ,q 

145.1 

4406.2 

Aro.O 

3 HO , 0 


1^7?..= 

P-f 52.-^ 

154.2 

3683 .6 

: .-. 

3H : , 3 


l-^'i. 3 

2rlFC.7 

1 ‘■ 3 . 3 

3C-:6. 4 

A J . 

3 *5 'v . C 


13 l=.4 

28‘’27.7 

172.4 

2 4 A 3 . A 

ACO.O 

HOO. 0 


1062.3 

2:491.8 

151.4 

1945.9 

'‘Cv..: 

H23.C 


812.6 

26167.2 

.90.5 

xl? 4 . 9 

^ • V# 

H40 .0 


^9^ , 4 

24856.6 

199, A 

110«. 2 


A-3J.0 


4^'^.= 

2-669.0 

20?. 7 

781. c 

CT^CULi' PfiYLf^AC 

= H6C.3 

LRS» 

209.0 '<C- 

cifculap v=l 

.» 24540. e 

F^S 

E5CADC PAifLC?AO 

a ■’9." 

LBS« 

36.1 >'G 

9SCAPE VEL. 

a 34718.6 

PPS 


Hi • >. 

3 =C.( 
2 2 C .L 
30C .C 

Av : ,c 

'• 00 . C 


<^1 .,;-iT C 4TI 

" , 'C-k A^ EL 0 ° E 
CA3= =AGc ^ 


- A7 - 



:v-Li:u u^'ITS 


71 IT ( • M mt a ^ j r i_ c 

LL«PTi:‘.L rv I'' 


VrtG^f C J P A T X^^> 

cwrc;^-^ Ef'P 


'* A S “ 


1 


o jr-t 


^E^IC-EE 4LT.S 


s CPA 


k y 


j. t:50O|_ ^ XT'"' ''ATA (CHECK FT= aCCL."^ACY) 


L <■ '. j 

T'' JECYT'^h 


V r L . . = = S 


- - ^ . c, 

w. ^ 1 

3 ’ \ 7a . *3 

1 = ' . ; 

■? r A ? 7 . 7 

:i c ^ \ 

2dar'.l 

30-3. •. 

’CA^.A.Q 

Atc .C 

?AA34. a 

50^.0 

? 29 5A , 7 

tij,Z 

ei'763.5 


P=3 IGcE 

P 4YLQA0 

APGGE f 

P I G £ c 

P AYLCAO 

APOGEE 

ALT.pN.YI. 

U T . , L 3 S 

ALT. .N.'^I . 

VEL.pFOS 

jt.,kg 

ALT.»KM 

EOC.O 

EO.C 

1AA67E. 2 

32 82 P .2 

3^.3 

257938.4 

? •• C . r- 

1 00 . 0 

282A7.9 

327fei.l 

45 .4 

7J835.1 

f j ■ j . - 

IPO.C 

212 AS .1 

31827.5 

54.4 

30242.3 

? 3 0 * 0 

1 c 

14311.1 

21 wCA ,6 

63.5 

265C4.1 

5C0.C 

160.0 

1G482. d 

3u2E7.9 

72.6 

19413.7 

f 0(3.C 

180. C 

8033.6 

29369.7 

81.6 

14878.3 

‘=C0.C 

2 C 0.0 

^^ 54,9 

289A2.5 

9\j • 7 

11750.2 

? V ^ • -I 

2 20.0 

5.3U.6 

28367.1 

QQ.- 

0501.0 

0 . ^ 

2 ^ 0.0 


3 ^ ? 3 • ? 

1 J F . c 

7725.3 

c - 

^ * V • V 

2^0. C 

3 ^ ^ • 7 

47^18.5 

117.9 

. 5 

f c 0 • c 

2 8 w . 0 

28 i6. A 

3 1 B ^ tj • ^ 

127.0 

521^.0 

E JC.C 

C'J*J • 

2212. d 

2 A 4 4 , Q 

*26.1 

42«A.e 

r j ■/ 1 

3 2^.0 

I s 4 2 . p 

3 

*45 .1 

3 A P 6 . 0 


3^3,0 

1 i ^ ^ 

t • ^ 

154.2 

P'^Or.O 

r — » ♦ ', 

?^J.O 

X • 2 

Z: 3.. :.3 

l‘-3.2 

2215 .A 


3 S j . 0 

^ . X 

33.C 

172. A 

17Ci .2 

C * '> 

J • • 

A ' n 

# w 

::7Z.f 

74^94 , 4 

1^1.4 

12A5 .A 

Cl=>CI'LAS P4YL0A 

A15 . A 

L^S= 

CI^CULA^ V-L 

.= 2^234,5 

p P 5 

F3rA°t D^VL‘1'‘0 

=« 72.2 

LBS* 32.7 K9 

ESCAPE VEL. 

= ’ 4775 ,*- 

FP8 


- A8 


U1 H . Ui X ^ f-* rr> r~l HJ C\j O in til o 



SC"UT V“HICL = 

£LLl®TICtL C.?'=5IT 3 = PF'I=>MNC= c^GLrSH IJMT$ 


VQUGWT CO^pOPATIC'^ 
P90G5A’’' FLI’P? 


CAS? 


PAG? 


PPPiC-c? ALT.= ^CC»C N.MI, »llll.Z 
r TE=P'-L AT OATA (C’rIcCK pr-B ACCURACY) 




■*T.,L = S 

Vf=L .^PPS 


3f 5 = 1. A 


hG.O 

i ?C ,•< 

2=716.2 

22 

27A2P .5 

300 .0 

25A28.9 

^uC.C 

23^69.2 

?nc.c 

2138A.1 

^ 0 0 # c 

2'J'5A7,C 


PPPIGEE 

PAYLOAD 

APOGEE 

PERIGEE 

PAYLOAD 

APOGEE 

ALT.,v|.M. 

WT.»LBS 

ALT.>N.MI. 

VEL.»FPS 

4T.,KG 

ALT.,Kf< 

600.0 

90.0 

69326.5 

329A6.2 

36,3 

128396,4 

^GC.O 

l&O.O 

28136.5 

■ 31870,1 

95.4 

52108,8 

600.0 

T20.0 

16998.6 

30928 .1 

54.4 

31481.4 

tcc.o 

1^0.0 

11323.6 

30097.1 

63.5 

21897,3 

60C.C 

160.0 

8782.3 

29342.6 

72.6 

16264.0 

□ jG.O 

180.0 

6763. C 

26646,8 

81.6 

12525.1 

6C0.0 

200.0 

5343,2 

28C11.8 

= 0.7 

9895.7 

ACQ.O 

2 20.0 

A288.7 

27A2P.5 

99. S 

7942.6 

6 90.C 

260.0 

-- 345A,fc 

26676.3 

108.9 

6398.0 

60C.0 

260.0 

2789.7 

26361.7 

117.9 

5166.4 

6C0.C 

28j.O 

2246, c 

256P0.6 

127.0 

4161.3 

600. C 

30C.0 

1794,9 

2542a. o 

136.1 

3324.2 

r^O.C 

320.0 

1405.5 

249=5.4 

145.1 

26C2.9 

6 C V . C 

340.0 

1-J70.O 

2^55 4. 

15A.2 

1983.3 

600. C 

360.0 

7fic. ? 

241=9.4 

163.3 

1^45.2 

CI°CIJLAC PAYLOAD 

= 374.0 

L6S= 169.7 Kr 

tTPCULAp VrL 

.= 23035 .C 

FPS 

ESCAPE PAYLOAD 

= 64.9 

L3S> 29.4 KG 

ESCAPE VEL. 

» 338A9.2 

FPS 


- A9 - 



S''?UT LiUNCM VPHiCLr 

f*LL:DTICAL PEPPnp -lAt^ce t^';LISH 'JMITS 


VCUGHT CO??OPATI^^J 
DCCG9AM ELOP‘D 


CASE 


D iCi 


PE=TC-cE ALT.* 7C0.C N.yi. ai?Qf..6 


l: TEPP?L AT I'^^ DATA (CHECK ACCURACY) 


ylt-ad 

TVJFCTIf^\ 

. >L'^S 

V.-L. •=p$ 

^ V t 

P^^lQ.a 


3iAO^, 5 

15S.L 

2a?Ai.4 

E 2 C • C 

26469.3 

3-0. C 

24426.7 

40C.C 

22416.0 

= 03.0 

2C771.8 

fOC.C 

19368. ? 


PERIGEE 

PAYLOAD 

APOGEE 

REP IGEE 

PAYLOAD 

APOGEE 

ALT.»N.«I. 

WT.,L«S 

ALT./U.f-I. 

VEL..FPS 

WT.#KG 

ALT.»KM 

700. C 

60.0 

489735.2 

33208.7 

27.2 

906989.6 

700. C 

80. G 

43635.5 

32057.2 

36,3 

80812.9 

700. C 

100.0 

21579.5 

30972 .C 

45.4 

39065.3 

70C.0 

120.0 

13774.0 

3CC20.3 

54.4 

25509.4 

700.0 

140.0 

9786.5 

29170.4 

63.5 

19124.5 

700.0 - 

160. C 

7323.8 

28414.7 

72.0 

13563.7 

7C0.G 

180.0 

5F-37.1 

27708.6 

81.6 

10439.9 

TOl.C 

2C0.0 

4425.0 

27063.2 

9C. 7 

8195.1 

700.0 

220. C 

351C. 3 

26469.3 

90.8 

6501.0 

7C0.C 

24j.0 

2778.2 

250c6. 4 

103.0 

5145 ,2 

700.0 

260 .C 

215B.8 

25361.0 

117.0 

4053.7 

700. C 

280. C 

17v,4.C 

24869.1 

127.1. 

3155.8 

70C.0 

300.0 

1297,7 

2^426.7 

136.1 

24C3. 3 

7f 0.0 

3 20. C 

Q4r . 2 

23983 . 1 

145,1 

1752.4 

CIRCULAR PAYLOAn 

» 336.2 

LPS* x52.5 XC 

CIRCULAR VEL 

.= 23644. 

4 PBS 

FSCAPE PAYLQAP’ 

» 57.7 

L8S» 26.2 KG 

ESCAPE VEL. 

» 33438 . 

3 F PS 


- AlO - 



SC3UT LAU*^'CH 7?HICLE VnUGHT CTRorRATICN 

ELLIPTICAL CRTIT p<=R FORJ^AHCS ENGLISH UNITS pRTGRAN CLCPE 

C6>c 1 PAGF sj 


P55IC-SE ALT.= POC.O » 1 A 6 1 . t- _ 


INT = '- POL iTIf‘'^ DATA (CHFC< F«^P ACCUPACY) 


paYLHAP 

ir jpcTioN 

■-■T. .L^S 

V?L.,=0S 

P<T 

3 ?® 3"^. 9 



1 3 */ • v» 

27S47. 3 

22i,i 

35AP3.8 

3C*C .C 

23391.9 

ACC.C 

21309,7 

. 501. C 

19539.7 

fcOC.C 

18103.9 


ofrigEF 

PAYLOAD 

APOGEE 

PERIGEE 

PAYLOAD 

APOGEE 

ALT.,N.NI. 

WT.tLPS 

■ ALT.,N.:^I. 

VEL.fFPS 

WT.#KG 

- ALT.#KM 

5CC.C 

60. C 

103653. C 

324C6.1 

27.2 

191965.3 

flOC.C 

80.0 

30521.5 

31153.0 

36.3 

56525.9 

80C.C 

ICO.O 

16903.2 

30055.9 

45.4 

31304.8 

eoc.c 

120. 0 

11191. fe 

29C92.1 

54.4 

20726.8 

soo.o 

140,0 

OU54.9 

28239.2 

63.5 

14917.7 

POO.O 

160.0 

6038,2 

27462.9 

72.6 

11182.7 

800.0 

130,0 

4820.3 

26745.5 

81 .6 

8556.8 

ecf .c 

200. 0 

3531,8 

26088.9 

90.7 

6633,4 

8C0.C 

220,0 

2736.7 

25483.8 

99.8 

5160.9 

ecG.c 

2 40 . C 

2145. u 

24QC9.0 

ioe.q 

3968,8 

800.0 

260, C 

Ifl^.c 

24371,5 

117. G 

3000.0 

flOC.C 

230. C 

1166,2 

23867,2 

127.0 

2196.8 

fcp.c 

3C0.0 

«2w.2 

233^1,9 

136.1 

1519.0 


cifculap payload* 

301.2 

L8S* 13o.o 

KG 

CIRCULAR VEL.» 233^4.2 

FPS 

FSCAPE PAYLGAH s 

50.6 

LBS* 23.0 

KG 

ESCAPh V5L. » 33042,0 

FPS 


- All - 
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‘^f'CUT L^L'NCi^ VEHICLf 

PLLI^TICAL rO'^.IT PEPFn=>'^ANC= ENGLISH UNITS 


VCUGHT COPPGPATICN 
PP^cpA^1 plOPE 


GAS. 


P^r-: 


PFPIGE': ALT.= QQC.O N.MI. =lc't6.? •<■/> 


1‘JTEPpOL ATI'^N CATA (CHECK F0p ACCURACY) 

payL'^AD I'JrCTIC''' 

UT.*L=S VEL..=PS 


A " , 

Irf .C 
220. C 
2 0C.C 
AO'J.C 
' 5C0.C 

1^00 . L 


32335. A 
2Q551.5 
Pbb74.3 
2R A 57 , 5 
223CJ. A 
2C120.2 
13315.6 
16730.8 


PEPIGEE 

PAYLOAD 


APOGEE 

PERIGEE 

paYload 

APOGEE 

ALT.fN.MI. 

WT .f LBS 


ALT.f N.MI. 

VEL.f EPS 

W T . » K G 

ALT., KM 

«500.0 

60.0 


53990.7 

31A90.3 

27.2 

99990.9 

300. C 

8 0.0 


22A7G.3 

3022A.4, 

36.3 

A1616.1 

900.0 

100.0 


13372.1 

2911A.8 

A5.A 

2A765.1 

90C.C 

120.0 


9069.3 

2S138.5 

5A. A 

16796.3 

90C.0 

lAO.C 


656A.5 

27272.8 

63.5 

12157. A 

900.0 

160. G 


A395.8 

26A82.7 

72.6 

9067.1 

°CC.O 

130.0 


3693.8 

2575C.5 

81.6 

68AC.9 

900.0 

2C0.0 


2797.7 

25C7S.5 

°C.7 

5181. A 

3tC.D 

220.0 


2102.6 

2AA57.5 

Q9 .8 

389A.1 

900. C 

2^0.0 


153A.7 

238<^6.9 

108.9 

28A2.2 

900. C 

260. C 


1^69.2 

22313.2 

117 .c 

1980. 1 

CIRCULAR PAYLQAO 

= 268 .5 

LBS = 

121.3 KG 

CIRCULAR VEL 

.a 23003.7 

pps 

ESCAPE PAYLOAC 

a A3. 5 

LBS« 

19.7 KG 

ESCAPE VEl. 

a 32659.5 

F°S 
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SC'.I-'T L^I)^JCH Ve-fTCLF 

-'LL^PTICiL C'PPIT 0 = -: Pn<5 MfinC 5 !:^GLISH U^'ITS 


vruG-'T GG 

?prG 


CASS 


PPPIGPS ALT.»1C')0.G N.mi, 


* 1 5 5 P , G 


:^ T‘=p?GL aTTG'“ TATt (CHFCK Pnp ACCURACY) 


0 A V L A f' 

injection 

wT. .LPS 

VEL. ,EPS 


3 2 : 7 . 4 

. L 

2 = r ^3 .2 

1 f • c 

? 5 « 33 . c 

2ZZ • ^ 

Z3371.C 

300# C 

? : 1 3 4 . 1 

^CG.C 

16353.2 

?CC,0 

16927.3 

ejc.c 

15225.5 


pPpIGbP 

PAYLOAD 

APOGEE 

PERIGEE 

payload 

ALT.,N .MI. 

WT.,LBS 

ALT.,N.mi. 

VEL ., EPS 

JT.,i^G 

1000.0 

^0.0 

247190. C 

320^,7.4 

18.1 

luOO.O 

60. C 

34408.4 

30546.5 

27.2 

iCOO.C 

60.0 

169C1.5 

2«264,4 

36.3 

lOCC.O 

ICO.C 

1U580.3 

28136.1 

45.4 

ICOO.C 

120. C 

7263.8 

27144.7 

5 4.4 

ICOO.C 

140. 0 

5238.8 

26261.4 

63.5 

ICOO.C 

1 6 0 . U 

3849. C 

25453.2 

72.6 

1 0 0 u .0 

180. C 

282 7. C 

24702 .4 

81.6 

icoc.c 

203. 0 

2053.1 

24C11.3 

= 0.7 

IOC 0.0 

220.0 

1445 . 5 

23373 .0 

9=.8 

CIPCULAP pAYLOad 

» 237.7 

LBS» 137.6 '<C- 

CIRCULAR VEL 

.» 22632 

ESCAPE PAYLOAH 

= 36.5 

L3S= 16.6 KG 

ESCAPE VEL. 

= 3226= 


sanotTir^ 
ELOPS 
PAGE IG 


APOGEE 
ALT., KM 


',5 77 95.8 
6372A.3 
31A68.9 
1959^.8 
13452.5 

9702.3 

7128.4 

5235.5 
3302. 3 
£C77.1 


pPS 


9 PPS 
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APPENDIX B 

FORTRAN CODE LISTINGS 


DR nr, RAM EL0PE( INRUT, OUTPUT# Tapes- INPUT»T A PE6 -OUTPUT# 

TAPE1#TAPE2#TAPE3» TAPES) 

PROGRAM ELOPE CONVERTS PERIGEE ALTI TUD E-VELOC I TY-PAYL OAO DATA 
TO PERIGEE ALTITUDE-APOGEE ALTITUDE-PAYLOAD DATA PARAMETRICALLY 


*t^¥************^^**>l^* INPUT DATA 


ALT^AX 

ALTMIN 

APOGEF 

DELALT 

DELWT 

FACT 

lOPT 


IPLOT 

IPRNT 

IRAO 

KASE 

LPLOT 

PEPIGE 

PLABELl 

-PLABLE9 

PTITLEl 

-PTITLEA 

PEARTH 

R1-P15 


MAXIMUM PERIGEE ALTITUDE USED FOR PARAMETRIC 
OUTPUT, UNITS ACCORDING TO lOPT. 

MINIMUM PERIGEE ALTITUDE USED FOR PARAMETRIC 
OUTPUT. UNITS ACCORDING TO lOPT. 

APOGEE ALTITUDE OF SINGLE ORBIT. UNITS ACCORDING 
TO lOPT. INPUT WHEN IQPT-3 OR A. 

INCREMENTAL PERIGEE ALTITUDE USED FOR PARAMETRIC 
OUTPUT. UNITS ACCORDING TO lOPT. MAXIMUM NUMBER OF 
ALTITUDE POINTS IS SO. 

INCREMENTAL PAYLOAD USED FOR PARAMETRIC OUTPUT. 

UNITS ACCORDING TO lOPT. 

FACTOR TO INCREASE OR DECREASE SIZE OF LOG PLOT. 

PLOT IS PLACED ON FILE PLT3. 

DATA OPTION 

-1 INPUT ALTMAX, ALTMIN, DELALT# DELWT# WTHAX, AND 
WTMIN IN N. MI.# LBS.# AND COMPUTE PARAMETRIC 
DATA IN N. MI. AND LBS. (1 BUILT-IN) 

-2 INPUT ALTMAX# ALTMIN# DELALT# DELWT# WTHAX# AND 
WTMIN IN KM AND KG AND COMPUTE PARAMETRIC DATA 
IN KM AND KG. 

-3 INPUT APOGEE AND PERIGEE IN N. Ml. AND COMPUTE 
SINGLE ORBIT ONLY. 

-A INPUT APOGEE AND PERIGEE IN KM AND COMPUTE 
SINGLE ORBIT ONLY. 

A NON ZERO VALUE PRODUCES A PLOT OF ALTITUDE VS 
VELOCITY. LOCAL PLOT FILE NAME IS PLT2. (0 BUILT-IN) 
FREQUENCY OF OUTPUT CONTROL. DATA IS COMPUTED AT 
DELWT INTERVALS FROM WTMIN AND PRINTED EVERY IPRNT 
DATA POINTS CALCULATED. (1 BUILT-IN) 

DATA OPTION 

=1 INPUT P1-R15 AS RADIUS IN FEET (1 BUILT-IN) 

-2 INPUT R1-R15 AS RADIUS IN N. MI, 

-3 INPUT R1-R15 AS ALTITUDE IN N. Hj, 

CASE number 

NON-ZERO VALUE PRODUCES A SEMI-LOG PLOT OF APOGEE 
AND PERIGEE ALTITUDE AS A FUNCTION OF PAYLOAD WEIGHT. 
LOCAL PLOT FILE NAME IS PLOT. (0 BUILT-IN) 

PERIGEE ALTITUDE OF SINGLE ORBIT. UNITS ACCORDING 
TO lOPT. INPUT WHEN IOPT-3 OR A. 

LABELS PLACED ON UPPER RIGHT SIDE OF APOGEE-PERIGEE 
PLOT. 30 TITLE CHARACTERS EACH. 

TITLES PLACED AT TOP OF SEMI-LOG PLOTS. 

AO TITLE CHARACTERS EACH. 

EARTH RADIUS USED TO COMPUTE ALTITUDES# FEET. 
(209257A1. BUILT-IN) 

TABLE OF RADIUS OR ALTITUDE (ACCORDING TO IRAD) FOR 
EACH WEIGHT. INPUT IN INCREASING ORDER. MINIMUM OF 
A AND MAXIMUM OF 10 VALUES PER TABLE. ENTER 0 AFTER 
LAST VALUE IF LESS THAN 10 VALUES ARE INPUT. 

MINIMUM OF A TABLES. 
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TITLE! « TITLE POINTED AT THE TOP OF EACH PAGE, ?2 TITLE 
CHAOACTEPS. 

TTTL52 » TITLE POINTED AT THE TOP OF EACH PAGE. 72 TITLE 

characters . 

V1-V15 « table of inertial VSLOCITY fop each WEIGHT. INPUT IN 

ORDER OF R1-R15. MAXIMUM OF 10 VALUES PGR TABLE. 
WEIGHT • TABLE OF PAYLOAD WEIGHTS/ IN LBS. CORRESPONDING RADIUS 
AND VELOCITY MUST BE INPUT FOR EACH WEIGHT. MINIMUM 
OF A VALUES AND MAXIMUM OF 15 VALUES. 

WTMAjr . MAXIMUM RAYLOAD WEIGHT USED FDR PARAMETRIC OUTPUT. 
UNITS ACCORDING TO lOPT. 

WTMIN > MINIMUM PAYLOAD WFIGHT USED FOR PARAMETRIC OUTPUT, 
UNITS ACCORDING TO lOPT. 

XINC » INCREMENTAL VALUE OF X AXIS MAJOR DIVISIONS/ LBS 
OR KG/ ON SEMI-LOG PLOT. THERE ARE SEVEN MAJOR 
DIVISIONS ON THE X AXIS. (100. BUILT-IN) 


implicit REAL(A-H/0-Z) 

COMMON /BLK1/TITLE1(8)/TITLE2 (6)/TITLE(l6)/LABEL(36) 

DIMENSION NC (15 ) / WW ( 15/ 50 ) / W6 16HT ( 15 ) 

dimension 9{10/15)/R1(10)/R2(10)/R3(10)/R4(10)/R5(10)/ 

A R6(10)/R7(10)/R8(10)/R9(10)/R10{10)/R11(10) / 

A R12(1C)/R13(10)/R14(10)/R15(10)/VI(10)/WI(10)/ 

A V(10/15)/ Vl( 10)/V2(10)/V3(10)/VA(10)/V5(10)/ 

A V6(10)/V7(10)/V8(10)/V9(10)/V10(10)/V11( 10) / 

A V12(10)/V13(10)/V1A(10)/ V15(10)/ 

A H(10/15)/HH(50)/VV(15/50) 

EQUIVALENCE ( Rl(l)/R(l/ !))/( R2(l)/R(l/ 2))/( R3(l)/R(l/ 3))/ 

A ( R4(l)/R(l/ 4))/( R5(l)/R(l/ 5))/( R6(l)/P(l/ 6))/ 

A ( R7(l)/R(l/ 7))/( R8(l)/R(l/ 0))/( R9(l)/R(l/ 9))/ 

A (R10(l)/R(l/10) )/ (Rll(l)/R(l/ll) )/ (R12(l)/R(l/12) )/ 

A (R13( 1)7R(1/13) )/ (R1A(1)/R(1,1A) )/ (R15(l) /R(l/15) )/ 

A ( Vl(l)/V(l/ !))/( V2(l)/V(l/ 2))/( V3(l)/V(l/ 3))/ 

A ( VA(1)/V(1/ ^))/( V5(l)/V(l/ 5))/( V6(l)/V(l/ 6))/ 

A ( V7(l)/V(l/ 7))/( V8(l)/V(l/ 8))/( V9(l)/V(l/ 9))/ 

A (V10(l)/V(l/10) )/(Vll(l)/V(l/ll) )/(V12(l)/V(l/12) )/ 

A (V13(1)»V(1/13))/(V1A(1)/V(1/1A))/(V15(1)/V(1/15)) 

DATA TITLE1/TITLE2/16*10H / 

DATA TITLE. LABEL/52+10H / 

data FTNM/6C76.11549//GM/l,A076576E16//XKG/,45359// XKM/1.852/ 
DATA NDIM/10/ 

namelist /INPUTD/IRAOz IOPT/ WEIGHT/ ALTMIN/ ALT MAX/ DEL ALT/ 

A WTMIN/ WTM AX/ DEL WT/R EARTH/ APOGEE/ PER I GE / R 1/ R 2/ 

A R3/RA/P5/R6/R7/R8/R9/R10/P11/R12/R13/R1A, R15/ 

A V1/V2/ V3/V4/ V5/V6/V7/ V8/ V9.V10/V11/V12/ V13/ 

A V1A/V15/LPLQT/IPRNT/ IPLOT/XINC/KASE/FACT 

C 
C 

C INITIALIZE DEFAULTS 

APOGEE-0. 

FACT»1. 

IOPT-1 

IPLOT-0 

IPRNT-1 
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I«A0«l 

K&SE-l 

LPLf^T«0 

PPPIG?»0. 

PEAPTH»2C9257^1. 

XINC-13C* 

C 

C PEAO INPUT DATA 

Ij CONTINUE 
CALL INPUT 
PEA0( 3>INPUT0) 

UNIT 1 USED FOP SEni-LOG PLOT 
UNIT 2 USED FOR SE'"I-L0G PLOT 
UNIT 3 USED FOR LINEAR PLOT 
REWIND 1 
RSyiNr 2 
REWIND 3 

deternine number of payload weights 
DO 20 I-l»15 

IF (W£IGHT( I) .EQ.O.) GOTO 30 
NWT»I 
20 CONTINUE 
30 CONTINUE 

DETERMINE NUMBER OF DATA POINTS PEP PAYLOAD 
00 60 I»1,NWT 

00 40 J«1>N0IM 

IF (R ( J> I ) .EQ.O. ) GOTO 50 

40 CONTINUE 

NC(I)»J 
GOTO 60 

50 NC(I)«J-1 

60 CONTINUE 
IPAGE-1 

CONVERT APPROPIATE INPUT TO ENGLISH IF METRIC 
GOTO (70»80,90.100), lOPT 
70 CONTINUE 

HMIN-ALTMIN 

mmax-altmax 
DH-OELALT 
WMIN«WTMIN 
WMAX-WTMAX 
0W»DELWT 
HAMAX-1.E6 
GOTO 110 
80 CONTINUE 

HMIN-ALTMIN/XKM 
HMAX-ALTMAX/XKM 
DH-OELALT/XKM 

wmin»wtmin/xi^g 
wmax-wtmax/xkg 

DW»OELWT/XKG 
HAMAX-1.E6/XKM 
GOTO 110 
QO CONTINUE 
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Ma. apogee 

HO «OE'? IGE 
GOTO lie 
*C/j comtimue 

HAaAPOGE6/XKM 
HP-PER IGE/XKM 
110 CONTIMIJE 

5E-PEARTH/FTNM 

CO’^VERT INPUT ALTITUCES OR RADII TO ALTITUDE# NM 
GOTO (12C#150# 16w )# IRAD 
1?C CHNTINUE 

DO lAC I»1»NWT 

DO 13U J«1#N0IM 

H(J,I)=R(J,I)/CTNM-R6 
IF (5 (J#I) ,EQ.C. ) H(J#I)-0. 

130 CONTINUE 

lAO CONTINUE 
GOTO 210 
150 CONTINUE 

00 170 I*1#NWT 

DO 160 J«1»NDIM 

H(J#I)»R(J»I)-RE 

IF (R ( J> I ) .EO.O. ) H(J#I)-0. 

160 CONTINUE 

170 CONTINUE 

GOTO 210 
180 CONTINUE 

DO 200 I-1»NWT 

DO 190 J-1#NDIN 
H( J»I)*R(J#I) 

190 CONTINUE 

200 CONTINUE 

210 CONTINUE 

00 230 I«1#NWT 

IF (I.EQ.l .OR, I.EQ.7 .OR. I.EQ.13) WRITE (6»A00) 
WRITE (6,420) WEIGHT(I) 

WRITE (6,A10) (H(J,I),J-1,N0IM) 

WRITE (6,430) 

WRITE (5,44u) (V(J,I),J-1,NDI'1) 

DO 220 J*1,NDIM 

WRITE (3) H(J,I),V(J,I) 

220 CONTINUE 
230 CONTINUE 

IF (IPLQT.NE.O) call PLOTLR (TITLE1,TITLE2,REARTH) 

IP (I0PT.GE.3) GOTO 370 
C 

C CALCULATE NUMBER OF PERIGEE ALTITUDES 

NH-INTI (HMAX-HMIN)/DH+1.5) 

NM-MINO (NH,50) 

C 

C CALCULATE VELOCITY FOR EACH PERIGEE AND PAYLOAD LINE 

HNEW-HMIN 
DO 25C x-1,NH 
MH (K ) »HNEW 
no 240 I»1,NWT 
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CALL INTEP(HNEW,VEL#NC(I),V(l>I)fH(l,I)) 
VV{ I,K )-VEL 
WW(I,K )-WEIGHT( I) 

9A0 Ca^JTI^UF 

HNEW»HNEw+DH 
E50 CONTINUE 
C 

C CALCULATE ORBIT DATA AT EACH PERIGEE 

OG 330 K»i»NH 
NREC»C 

HHNET»HH(K)*XKM 

WRITE (6,A50) TITLE1»TITLE2»KASE» ip age 
NBITE HM(K)»mmmeT 

WRITE (6>A70) 

WOITE (6»AB0) (WW(I,K),V'/(I,K),I«1,NWT) 

V»ITE (6»AQ0) 

HF»HH(K) 

RF«(HF+RE)*FTNM 

C 

C REVERSE ORDER OF VELOCITY AND WEIGHT ARRAYS 

J»NWT 

no 260 I-1>NWT 
VI(I)»VV(J>K) 

WI(I)-WW(J>K) 

J«J-1 

260 CONTINUE 
C 

C CALCULATE PAYLOAD AT 1E6 APOGEE ALTITUDE 

A-{ (HAMAX + RE)*FTNM+RF )/2. 

VMAX-SQRT (GM*(2./RF-1./A) ) 

CALL INTER(VMAX>PL>NWTi WI»VI) 

IF (PL. LT. WHIN) GOTO 270 
NREC-NREC+1 
WRITE (1) PL^HAMAX 
270 CONTINUE 
C 

C calculate NUHBER of payload WEIGHTS 

NNW«INT( (WMAX-UMIN)/DW+1,5) 

C 

C CALCULATE ORBIT DATA AT PAYLOAD INCREMENT 

WNEW-WM IN 
DO 300 I*1»NNW 

CALL INTEP{WNEW,VEL,NWT,VV(1»K)#WW(1,K) ) 

IF (VEL.GT.VMAX) GOTO 290 

HA-(2. / (2. /RF-VEL**2/GM )-RF) /FTNM-RE 

IF (HA.LT.HF) GOTO 310 

IF (noD(I-1»IPRNT).NE.O) goto 280 

PLM»WNEW^XKG 

hanet«ha*xkn 

WRITE (6,500) HF,WNEW,HA,VEL,PLM,HAHET 
280 CONTINUE 

IF (LPLOT.EO.O) GOTO 290 
IF (HA.GT.HAMAX) GOTO 290 
WRITE (1) WNEW,HA 
NREC-NREC+1 
290 CONTINUE 

WNEW-WNEW+DW 
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^.00 CONTINUE 

C 

310 CONTINUE 

r CALCULATE CIRCULAi? AND ESCAPE ®AYLCAO AT PERIGEE ALTITUDE 

V?LC»SORT(G,N/PF) 

CALL INTER (VELC^WTCfNWTy WI, VI ) 

WTC^ET»WTC*XKG 

VELE-SORT (2. )*VELC 

CALL INTER(VELE»WTE»NWT»WI,VI) 

wTEiETaUTE*VKG 

WRITE (6»510) WTC>WTCMET,VELC*WTE»WTENET» VELE 
10 AGP«I°AGE+1 
IF (LPLOT.SQ.n; ) GOTO 330 
IP (WTC.GT.WYAX) GOTO 320 
■WRITE (1) WTC>HF 
NREC»NREC+1 
320 CONTINUE 
C 

WRITE (2) NREC 
330 CONTINUE 

IF (LPLOT.EQ.O) GOTO 39C 
C 

C CALCULATE CIRCULAR QPBIT LINE 

HP«HMIN 
M A«HMIN 
NREC-0 

340 CONTINUE 
J-NWT 

no 350 I»1>NWT 

CALL INTER(HP,VEL,NC(I)>V(1>I)>H(1>I)) 

VI ( J) »VEL 
WI(J)»WEIGHT(I) 

J«J-1 

350 CONTINUE 

RBAR»(riP+RE)+FTNM 
A«(2.^RE+HP+HA)/2.TFTNM 
VEL-SQRT (G'^* (2./RBAR-1./ A) ) 

CALL INTER (VEL^PL. NWT>WI, VI) 

IF (PL.LE.WMIN .OR. PL. GT. UMAX) GOTO 360 
NR EC-NPEC + 1 
WRITE (1) PL. HA 
360 CONTINUE 

ma-HA+DH/IC. 

HP»HA 

IF (HP.LE.HMAX) GOTO 34C 
WRITE (2) NREC 

CALL PLOTLG { IQPT. X INC. ALTN IN » DEL AL T, NH, 1 .. A 4. 4HPL0T ) 

CALL PLOTLG ( lOP T» X INC. ALTN IN, DEL AL T. NH, F ACT, 36. 4HPLT3 ) 

GOTO 390 
C 

C CALCULATE SINGLE ORBIT 

370 CONTINUE 
J-NWT 

DO 3B0 I-l.NWT - - 

CALL INTER (HP.VEL.NC ( I) ,V (1,1 ),H(1, I ) ) 

VI (J ) -VEL 

WI ( J ) -WEIGHT( I ) 
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J»J-1 
33J CONTINUE 

pa ^P« (HP+5E ) 

A» < 2.*‘?E+MO + HA) /2.*FTN.M 

\/EL-S'3RT (G'<* (2./P8AR-1./A) ) 

call inter (V6L»PL^NWT,WI, VI ) 

HPNET-HP+XKM 

MAMgT-HA^TKN 

OL*'«P L*XKG 

WPIT? {6/AEO TITLElf TITLi2#KASE> IP&GE 
WRITE (6>^6C) HP.HPMET 
-RITE (6»A90) 

'-"RITE (e»500) MP,PL»HA»VEL#PLN,HAW|£T 

39J CONTINUE 

K AS.= «KASS+1 

GOTO 10 


'•♦00 FORMAT 
^♦10 FORMAT 
420 POPMAT 
430 FORMAT 
440 FORMAT 
450 FORMAT 


4fcu FORMAT 
470 FORMAT 


480 FORMAT 


(IHl) 

( (15X,5(F10.3))) 

( / />15X,«WEIGHT(LBS. )«^,F8.1/15X, ♦ALTITUDE (N. MI, )*) 
( 15X,^VEL0CITY(FPS) + ) 

( (15X,5(F10.1) ) ) 

(♦!♦# 8X,8A10»T68»#V0U6HT COR PORAT lON^# / 
9X,8A10,T73, ♦PROGRAM ELOPE^>/ 

T68#^CASE *tlZ, 

T79»^PA6E ♦»I2//) 

( 25X#^P£RIGEE ALT.-^»F6.1>^ N.MI. ■♦jFO.l#^ KM^//) 
(15X, ♦INTERPOLATION DATA tCHECK FOR ACCURACY )♦,/ / 
13X,^PAYL0AD+>6X>^INJECTI0N^,/ 

13X »♦«!.# LB S^#7X#^VEL.#FPS+,/) 

(5X»2F15.1) 


400 FHRMAT ( ///13X,^PERIG£E^#9X,^PAYL0A0+#8X>+AP0GEE+#8X,+PERIGEE*> 
A 3X,4PAYLnA0^»4X,^AP0GEE*^ / 

A 13X,^ALT.,N.MI.^,6X,^WT.,LBS*»6X, ♦ALT.,N.MI,^, 

A 5X,^VEL.» FPS^»4X,^WT.»KG^#3X»*ALT.>KM^# / ) 

500 FORMAT ( 5X » 4F1 5 . 1 # 2F 10 . 1 ) 

510 format ( ///13X, ♦circular PAYLOAD-^, F6.1,+ LBS-^»F6,1,^ KG^» 

A SXf^CIRCULAR VEL.-^»F8.1»^ FPS^W/ 

A 13X, ♦ESCAPE PAYLOAD ■♦,F6.1,* L3S-^,F6.1,^ KG^, 

A 5X, ♦ESCAPE VEL. •♦»F8.1#^ F^S^) 

END 
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SUBPQUTINc TbP'JT 

TMIS SU8P0UTIN6 '?E40S MODIFIED NAMELIST FORMATTED DATA. 
IT READS A CARD DN UNIT 5»WRITES THE CARD ON UNIT 6> 
-RITES The card IN UNIT 8 (FIRST 72 CHARACTERS ONLY). 
THE TITLE CARDS AS DEFINED IN THE DATA STATEMENT BELOW 
ADE NOT WRITTEN ON UNIT fl BUT THE DATA IS PLACED IN 
THP appropriate arrays P0« TRANSFER BACK TO THE CALLING 
PROGRAM. THE TITLE CARDS MUST BEGIN IN COLUMN 2 WITH 
NO SPACES. THE CALLING PROGRAM MUST BLANK THE TITLE 
ARRAYS^ CALL IN^UT AND P EAC (8» INPUTD ) . NAMELIST DATA 
'*UST BEGIN WITH SIMPUTO AND END WITH SEND, BOTH 


'>==GTNNINf: IN COLUMN 2. 
T'*PIICIT INTEGER (4-7) 
Cr»»MON /BLWl/ TITLE(6P) 
dimension CAPC(8), LINtdS) 
DATA LINE/IGH TITLEl- ,10H 

TITLE2- 

,10H 

PTITLEl 

, i'JH 

PTITLE2 

^ ,1CH PTITLE3 

,10H 

PTITLE4 

,10H 

PLA8EL1 

,10H 

PLA8EL2 

^ ,10H PLA0EL3 

,10H 

PLABEL4 

,10H 

PLABEL5 

,10H 

PLABEL6 

^ ,10H DLABEL7 

,10H 

PLA8EL8 

,10H 

PLABEL9 

/ 


DATA BLANK/ICH 

/ 







rewind 8 

WRITE (6,70) 

10 CONTINUE 

READ (5,110 CARO 
IF (EOF(5).NE.O) STOP 
WRITE (6,90 CARD 

BLANK COLUMNS 9 AND 10 
ENCODE (10, 8C, WORD) CARD ( 1 ), BLANK 
DO 30 1-1,15 

IF (WORD.NE.LINE(I) ) GOTO 30 
CARD READ IS A TITLE CARO 
IF (I.EO.l) J»1 
IF (I.EQ.2) J-9 
IF (I.GE.3) GOTO 20 
ENC00E(72,60,TITLE(J) )CARD 
GOTO 10 

2vT CONTINUE 

j-17+4*(I-3) 

=NC0OE(3O,5O,TITL5(J) ) CARD(l), CARD(2), CARD(3), CAPD(4) 
GOTO 1C 

30 CONTINUE 

BLANK columns 73-80 OF DATA CARD 
ENCODE (1G,1CD,CARD(8) ) CAR0(8), BLANK 
WRITE (8,110) CARO 

IF (CARO(l) .NE.lOH SEND ) GOTO 10 

REWIND 8 
CALL OATE(DAT) 

CALL TIME(TIM) 

WRITE (6,40) DAT, TIM 
RETURN 
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cqcv*at 

( /////>ijX,*DATE 

A 


10X,*TIME 


FHRM4T 

( Rif AiOf AlOf A9 ) 

hj 

FC’^'^AT 

(52»7A1J) 

70 

format 

(IH*) 

-0 

CQRVAT 

( AP, A2) 

<50 

FP.RwaT 

( 10X,8A10) 

IOC 

format 

< A2f A9 ) 

110 

format 

EKD 

( 8A1C ) 


IS *, 6 . 9 / 
IS *,\ 9 ) 
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SUBRHUTINe INTER ( X » Y, ML'M, B » A ) 

SPCINO Ofir.ER INTERPOLATOR 

SELECT FGliP DATA POINTS CLOSEST TO X TO INTERPOLATE FOR Y, 
x»ImoE°ENOENT variable VALLE 

y-resulting deoenoent variable value 

L’^TaNO. OF ELF.NENTS IN A AND 8 
C BaARRAY OF DEPENDENT VARIABLES 

c A-ARRAY OF INDEPENDENT VARIABLES 

DI*^ENSION A(15)>B(15) 

C 

1*1 

IF (N'J‘^.EO.4) GOTO 30 

IF (MUM.LT.A) VRITE (6M0) NUW 

IF (NUN.LT.A) STOP 

C 

IF (X.LT.AO)) lal 

IF (X.GT,A(NUN-2) ) IaNUN-3 

IF (X.LT.AO) .OR. X.GT.A(NUM-2)) GOTO 30 

C 

LMTaNUH-2 
DO 10 K-AfLNT 

IF ( X.LT.A(K) ) GOTO 20 
10 CONTINUE 
20 CONTINUE 
I»K-2 
C 

30 CONTINUE 
XO»A( I) 

Xl-A( I+l) 

X2«A(I+2) 

X3aA< 1+3) 

Ylla( (X1-X)*8(I)-(X0-X)^P(I+1) )/(Xl-XO) 

Y21»( (X2-X)*B(I)-(X0-X)*B(I+2) )/(X2-XC) 

Y 31 a( (X3-X)*B( I)-(XO-X)*B (1+3) )/(X3-X0) 

Y22»( (X2-X)*Y11-(X1-X)*Y21)/ (X2-X1) 

Y32=( (X3-X)*Y11-(X1-X)+Y31) /( X3-X1) 

Y* ( (X3-X)+Y22-(X2-X )+Y32) /( X3-X2) 

R ETURN 
C 

AC FOR^«AT ( //iCX, ♦SUBROUTINE INTER - VALUES IN INTERPOLATION TABLE ■♦ 

13* MUST BE .G5. A+) 

END 


SIC 



o o -> o o o r> o o 


SI-BPOUTIME PLOTLC- I npT,XiNC,HMiN,OH,NUMH»FACT#ICAL/PPlLE) 
C T-iIS SUnPOUTINe G5MEPATSS A SEMI-LOG PLOT OF APOGEE 

C MTTTL'DE A5 A FUNCTION OF PERIGEE ALTITUDE AND 

C saYLHAO WEIGHT. 


INPUT DATA 

DM * PERIGEE altitude INCREMENT 

FACT « FACTOR FOR RELATIVE SIZE OF PLOT PRODUCED. 

hmim . mimmum perigee altitude 

ICAL « TWO DIGIT CALCCDE. 

lOPT * I IF ENGLISH UMTS 

■ 2 IF METRIC UNITS 
MUMM > number of perigee ALTITUDES 
PFILS « LOCAL FILS NAME OF PLOT FILE. 
i^IMC « A3CISSA MAJOR DIVISION INCREMENT 


C 

C 


C 

C 


IMPLICIT REAL(A-H,0-S#U-Z)f 
. INTE6ER(I-N,T) 

COMMON /BLK1/TITLE1(8)»TITLE2(8)>TITLE(16)#LABEL(36) 
DIMENSION X(lOOO)> Y(IOCO) 
dimension A(5)» LABELY{20) 


A 

A 

A 

A 


DATA LABELY/IOHIOO 
^ 10H600 

^ 10H2000 

1GH8000 
10H40000 
ICHIOOOOO 
10H600000 


#10H200 

>10H800 

>10HA000 

♦lOHlOOOO 

» 10H60000 

f 10H20G000 

>10H800000 


>10H400 
»10H1000 
f 10H6000 
>1CH20C00 
>10H80000 
«1CH400000 
/ 


data XKG/0. 45359/, XKM/1, 852/ 

ENCODE(5,5,CALCODE) ICAL 
F0RMAT(3HCAL,I2) 

CALL PLOTS (CALC0DE,0,PFILE) 
CALL FACTQR(FACT) 

REWIND 1 
REWIND 2 
NCYCLE-4 
A(1)»0. 

A {2)»AL0G10 (2. ) 

A (3)-AL0G10(^. ) 
A(4).ALQG10(6.) 
A(5)»ALOG10(e. ) 

YINT-0. 


DC 
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TICK M 

APKS 

ON 

Y AXIS 

AND ANNOTATE 

20 

J-1,NCYCLE 




DO 

10 I« 

1,5 





Y1-{FLQAT( 

J- 

1)+A( I) 

)42.5 


CALL 

PLOT 

(C 

.,Y1,3) 



CALL 

PLOT 

(- 

.05, Yl, 

2) 


<1— . 

1*(3. 

+ FLOAT( J ) 

) 


Y1*Y1 

-.05 





I3CD- 

LABELY ( 

5*J-5+I 

) 


NCHAR 

■2* J 





CALL 

SYMBOL 

(XI, Yl, 

.1, IBCD,C.,+NCHAR) 

CONTINUE 
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20 CGNTI^UE 

Y1-(FL0AT{NCYCLE)+A(1) )*2.5 
CALL °LQT (0.»Y1,3) 
call plot (-.05^Y1»2) 
vi»-,l^(3.+FLnAT(NCYCLE+l) ) 

Y1-Y1-.05 

CALL SY'IBOL ( X 1, Y 1> . 7H100GOCO» 0 . # +7 ) 


C 

C 


C 

C 


C 

C 


OPAV TICK -iACKS ON X AXIS AND ANNOTATE 
on 3J J«l,« 

Xl«FLr!AT( J-1) 

CALL PLOT (XI, G. >3) 

CALL °LGT (Xi,-.J5»2) 
xl=xi-,15 

FPN»X INT+ ( J-1) *X INC 

CALL NUH3ER(X1,-,2# .1»FPN,0.#-1) 

33 CONTINUE 


LABEL X AND Y AXIS 
IP lIOPT.EQ.i) 

A CALL SYMB0L(-.6A>2.A5».1>21HAP0GEE 
IP (IQPT.EQ.l) 

A CALL SYMB0L(2,3#-.35# .1>21HPAYLQAD 
IF (I0PT.EQ.2) 

A CALL SYMB0L(-.64,2.A5> ,1,18HAP0GEE 
IF (I0PT.EQ.2) 

^ CALL SYMBOL (2.3,-. 35» .1,24HPAYL0AD 


ALTITUDE-N.MI.,90.,+23) 
WEIGHT-PQUNDS,0. ,+23) 
ALTITUDE-KM,90., +20) 
WEIGHT-KIL0GRAHS,0.,+23) 


WRITE 4 LINES OF TITLE AND 9 LINES OF LABEL 
K «1 

Xl-7./2.-8.*.14 
Yl*10.-.14 
00 40 1-1,4 

CALL SYMBOL (X1,Y1, .14,TITLE(K),0.,30) 
Y1-Y1-.14-.035 
K-K + 4 
40 CONTINUE 


K-1 

X1-7.-.10+30.5 
Yl-8 . 5 
00 50 I-l,«3 

CALL SYMBOL ( XI, Yi, , 10, L AB EL ( K ) , 0 . , 30 ) 
Yl-Yl-.lC-.035 
K-K+4 
50 CONTINUE 


FILL ARRAYS FOR PLOTTING 
60 CONTINUE 

READ(2) NREC 

IF (E0F(2) .NE.O) GOTO 80 
DO 70 I-1,NREC 
REAO(l) WT,HA 


70 


IF (lOPT.EQ.l) 
IF (lOPT.EO.l) 
IF (I0PT.EQ.2) 
IF (I0PT.EQ.2) 
CONTINUE 


X( I)«WT 

Y ( I ) -AL0G10(HA )-2. 

X ( I ) -WT*XKG 

Y( I )»AL0G10{HA+XKM)-2. 
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DEFINE SCALE FACTORS 
X( M5 FC + 1 ) »XINT 
X(NREC+2 )«XINC 
Y (NREC+1) -C. 

Y(NPEC+2)»1./2.5 
C 

r ''RAW CURVE 

CALL LINE (X,Y,NREC#1>C,0) 

GOTT fo 

so CONTINUE 
C 

C WRITE additional LA^JELS IN RIGHT BQRDER QF ®L0T 

X1*F. 

Yl*6. 

HGT.,«-)35 
IF (lOPT.EO.l) 

A CALL SYMBOL ( XI, Yl,HGT#2^HPERtGEE ALTITUDE - N,MI .,0.,+24) 
IF (I0PT.EQ.2) 

A CALL symbol (X1,Y1,H6T,21HPERI6EE ALTITUDE - KI1»0.»+21) 

DO 00 I»1,NUMH 
Y1-Y1-HGT-.035 
FPNaHMIN + OH>*( I-l) 

CALL NUMBER(X1,Y1,H6T,FPN,0.,-1) 

00 CONTINUE 
Yl.Yl-,2 

CALL SYMBOL ( XI, Yl, HGT» 15HC IP CULAR ORB ITS,0. , +15 ) 

CALL PLOT (12. ,0., 999) 

RETURN 

END 
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SUBROUTINE PLTTLP (TITLE 1 #TITL 52 »R EARTH) 

THIS SUBROUTINE GENERATES A LINEAR PLOT OF PERIGEE ALTITUDE 

AS A fijnctton of velocity and payload weight, 

implicit R£AL(A-H, 0 -S,V-Z)> 

INTEGEPC I-N,T) 

DIMENSION TITLEK 8 ), TITLE 2 ( 8 ) 

DIMENSION X( 152 ), Y(i 52 ), Z( 150 ) 

DATA CN/ 2 . 5 A/,GM/i. 4076576 E 16 / 


CALL PLOTS( 5MCAL32»0,AHPLT2) 

REWIND 3 
SC2»SORT (2, ) 

C 

c ANNOTATE X-AXIS 

nchar=5 
no Iw J»1,2C 

Xl-FLOAT(J-l)*2./CM 
CALL PL0T(X1^0.»3) 

CALL PL0T(X1,-,Q5>2) 

X1-X1-.1+NCHAR/2. 

FPN»18000. + 10C>C,*( J-1) 

CALL NUMBER (XI,-. 2> ,1>FPN,0.»-1) 

10 CONTINUE 

CALL SYMBQL(5.,-.35, .1,23HINJECTI0N VELOCITY, FPS,0.,23) 
C 

C ANNOTATE Y-AXIS 

NCHAR-4 

on 20 J»l,13 

Y1-FL0AT(J-1)*2./CM 
CALL PL0T(0.,Y1,3) 

CALL PL0T(-.05,Y1,2) 

Y1-Y1-.05 

Xl»-.05-.l*NCHAR 

FPN-0.+iec.*(J-l) 

CALL NUMBER ( XI, Yl, .1,FPN,0.,-1) 

2J CONTINUE 

CALL SYM3QL(-.6,2.5,.1»22HINJECTI(3N ALTITUDE, NM,90.,22) 
C 

C WRITE TITLES 

X1-8.-.1A+A0. 

Y1-10.-.07 

CALL SYMB0L(X1,Y1,.1‘», TITLE1,0., 80 ) 

Y1-Y1-.21 

CALL SYM80L(X1,Y1,.14, T IT L £2, 0 . , 80 ) 

C 

C DRAW CIRCULAR AND ESCAPE VELOCITY LINES 

H»100. 

T-1 

3U CONTINUE 
X( I)»H 

Y( I)»SQRT(GM/(Xm + 6L76.115A9+REARTH) ) 

Z( I)»SQ2*Y(T ) 

IF (H.EQ.120C.) GOTO 40 

I-I + l 

H»H+20. 

GOTO 30 
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40 CONTINUE 
y ( I +i )»c. 

X ( 1+2) »10j./2.*C‘^ 

Y( T*1 )aieOOC. 

Y(I+2)»1000./2.*CM 

z(T+i)«ieooo. 

Z( I+2)»l0DC./2.*CN 
CALL LIMt( Y,y,I,l,0»0) 

CALL LINE(Z^X»I>1>0>C) 

PLOT DATA 
I»1 

tO CD'^TI^UE 

P£AD(3) Y(I),y (I) 

IF (EOFt 3) .NE.C) GOTO 60 
IF (Y(I).LT.O .OR. Y ( I ) .GT .1200. .OR. 

A X( I) .LT.18000. .OR. X ( I ) .GT .37COO. ) GOTO 5C 

I-I + l 
GOTO 50 
60 CONTINUE 
I-I-l 
Y{ I + 1)=«0. 

Y{I+2)-lC0./2.*CM 

X(I+1)«18000. 

X( I+2)»1C00./2.*CM 
CALL LTNE(X,Y,I,1,-1,1) 

CALL PL0T(17.,0.,999) 

RETURN 

END 
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APPENDIX C 


SCIENTIFIC DATA PROCESSING ROUTINE 
SUMMARY DOCUIffiNTATION 


IDENTIFICATION 

Title Elliptical Orbit Performance 

Routine No. 7031 Date Filed March 72 Security Class. U 

Responsible Engineer T. R. Myler 

Date Completed March 1972 Source FORTRAN Other 

Language; IV 

Key Words Orbit parameters, interpolation, C alComp plot 

RESOURCE REQUIREMENTS 

Typical CPU 5 sec Machine (s) CPC CYBER 175 No. Source Cards 710 

Core 60k (octal) Tape none Plot yes Graphics none 

DESCRIPTION 


Purpose ; To calculate elliptical orbit altitudes as a function of payload 
weight and generate a CALCOMP plot. 

Input ; Parametric data of insertion velocity as a function of insertion 
altitude and payload weight. Also selectors on desired plot. 

Output ; A table of apogee altitude and payload weight for each perigee 
altitude selected. Also the payload weight at the circular and escape 
condition. All data output in both English and metric units. CalComp plot. 

Functional Description ; Uses a 2-body orbit equation for orbit determination 
and a second order curve fit to interpolate for velocity and payload weight. 


DOCUMENTATION 

Vought Report 2-53030 /lR-5 2643 » "Elliptical Orbit Performance Computer 
Program" dated L June 1981. 
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